Sequential Generation of Matrix-Product States in Cavity QED 
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We study the sequential generation of entangled photonic and atomic multi-qubit states in the 
realm of cavity QED. We extend the work of C. Schon et al [Phys. Rev. Lett. 95, 110503 (2005)], 
where it was shown that all states generated in a sequential manner can be classified efficiently 
in terms of matrix-product states. In particular, we consider two scenarios: photonic multi-qubit 
states sequentially generated at the cavity output of a single-photon source and atomic multi-qubit 
states generated by their sequential interaction with the same cavity mode. 

PACS numbers: 03.67.Mn, 42.50.Dv, 42.50.Pq 



I. INTRODUCTION 

Triggered single-photon sources p], 0, Q have impor- 
tant applications in quantum communication and quan- 
tum computation and are subject of intense experimental 
investigation 0, 0, i, 0, H, i, EE El E3, El, El, EE El, E3 • 
Substantial progress is being made along various routes 
using single atoms 0,0.00, II, ions IH, molecules 
Till, quantum dots [fl2l. IT3L TT3 . EH or color centers 
id |l7[ . The basic operating principle common to these 
approaches is that a triggered single quantum emitter 
excites the mode of a cavity or a photonic crystal with 
a single-photon, which coherently leaks out into a well 
defined field mode. If the source is run in a cyclic fash- 
ion: sequentially initialized, loaded, and fired, the result 
is ideally a train of identical single-photon wave packets. 

Given the impressive experimental achievements along 
these lines, it is interesting to go one step further and 
address the following scenario. Assume that the source 
is not initialized after each step, but stays in some quan- 
tum state, which can in turn be correlated to the field 
state generated so far. What kind of multipartite quan- 
tum states can then in principle be created with such a 
sequential generation scheme? We answered this ques- 
tion in Ref. [18[ and proved that the class of sequentially 
generated states is exactly identical to a class of so-called 
matrix-product-states (MPS) 0, E El . These states 
play an important role in a completely different context, 
namely in the theory of one-dimensional spin chains [22[ , 
where they constitute the set of variational states over 
which Density Matrix Renormalization Group techniques 
are carried out [23|, 0, HBj]. This classification is sig- 
nificant in at least three respects: (i) It does not only 
apply to all the various types of single-photon sources 

but literally to any system where a multipartite quantum 
state is generated by sequential interaction with a source- 
or ancilla-system [26] . (ii) It provides a constructive pro- 
tocol for the generation of any desired quantum state. It 
allows one to decide whether or not a given state can be 
generated with a given setup and, if so, what sequential 



operations one has to apply, (iii) It stresses the impor- 
tance of the MPS formalism as it demonstrates that they 
naturally occur in vastly different physical contexts. 

In this work, we elaborate on the results of Ref. [l8[, 
providing more details and further applications. In sec- 
tion [III we fi rs t illustrate the basic idea of sequential gen- 
eration of quantum states by means of a cavity QED 
(CQED) single-photon source 0, H, Q. Then, we give 
a detailed derivation and extension of the main results 
in Ref. [18[ without reference to any specific setup. In 
section IIII( we show how to create important entangled 
multi-qubit photonic states encoded in photonic time 
bins and in polarization states. In section IIV[ we con- 
sider a microwave cavity interacting with a sequence of 
atoms flying through it, such as in Refs. [28|, |29[, gener- 
ating entangled multi qubit atomic states. 



II. SEQUENTIAL GENERATION OF 
ENTANGLED MULTI-QUBIT STATES 



In this section, we consider first the CQED single- 
photon source |4|, |5|, uj and show, with an elementary 
example, how the structure of MPS arises quite natu- 
rally in this setup. We then treat the sequential genera- 
tion scenario without referring to any particular physical 
system. First, we assume that an arbitrary source-qubit 
interaction is available in each step of the sequential gen- 
eration. In the second part, we restrict the interaction 
in a way which resembles the situation in current cavity 
QED setups for the generation of single-photon pulses. 
Finally, we discuss a more abstract scenario, a register of 
qubit s, to which nearest neighbor gates are applied se- 
quentially. This situation is again covered conveniently 
by the MPS formalism, even if the sequence of gates is 
applied repeatedly. 
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FIG. 1: A trapped D- level atom is coupled to a cavity qubit, 
determined by the energy eigenstates |0) and |1). After bipar- 
tite source-qubit operations, photonic time-bins are sequen- 
tially and coherently emitted at the cavity output, creating a 
desired entangled multi-qubit stream. 




FIG. 2: The hen and egg picture: The ancilla A couples 
sequentially to initially uncorrelated qubits B%. In the last 
step we require the ancilla to decouple from the entangled 
n- qubit state. 



A. Basic idea 

The CQED single-photon source utilizes the possibility 
to excite the cavity mode via an atom, which is trapped 
inside the cavity 0, H, Q . The photon is then coherently 
emitted through the cavity mirror. A photonic qubit may 
be defined either by different polarization states or by the 
absence and the presence of the photon. If we allow for 
specific operations inside the source before each photon 
emission, we will be able to create different multi-qubit 
states at the output (see Fig. [1]). 

Typically, a CQED single-photon source employs an 
effective two- level atom with hyperfine ground states \a) 
and I b) coupled via a Raman transition driven by the 
cavity mode and an external laser. The latter controls 
the population transfer 

|a, 0) -> cos(0i)|a, 0) + sin(0i)e^ 1 16, 1), (1) 

where |0) and |1) denote the absence and the presence 
of a photon in the cavity mode. This so-called time-bin 
qubit is then coherently emitted as depicted in Fig. dJ 

We define q = cos(^) and Si = sin(^)e*^ i for step i 
and repeat the procedure n times. We end up with 

\a) -> ci|a,0) + si|&,l) 

-> cic 2 |a, 0, 0) + 1 6) (ci5 2 |l, 0) + si |0, 1)) 

|6)[ci...c n _i|l,0,...,0) 

+ ci . . . c n - 2 s n -i |0, 1, 0, . . . , 0) + . . . 

+ ci5 2 |0,...,0,l,0) + 5i|0,...,0,l)], (2) 

where we chose c n = and s n = 1. Then, the atom 
decouples and the resulting photonic state is a W-type 
n-qubit state. This is just a particular example of a more 
general sequential generation scheme (see Fig. [2]). 

An ancillary system A (the atom) with Hilbert space 
Ha — & D (D=2) couples sequentially to initially un- 
correlated qubits Bi (the time-bin qubits) with Hilbert 
spaces TLb — C 2 . In every step we have a unitary time 
evolution of the joint system Ha ® Hb- Since each qubit 
is initially in the state |0), we disregard the qubits at the 
input and write the evolution in the form of an isometry 
V : Ha —> 7~(-A ® Hb- The photon-generating process 
discussed above is given by 

% = Ci\0)i\a)(a\ + 8i\l)i\b)(a\ + \0)i\b)(b\, (3) 



for step i and fulfills the isometry condition Vjjj Vj^j = 1 2 • 
The final state can then be written as 

|*> = V [n] ...V[ 1] | ¥ »j> ) (4) 

where \ipi) = \a) is the initial state of the atom. Since 
the atom decouples in the last step, we may also write 

I*) = \<PF)(<PF\V[n\---V ll] \<p I ) 

= \<Pf)\1>), (5) 

where \^>f) — \b) is the final state of the atom and 
the n-qubit photonic state. The goal in the following 
sections is to classify all achievable states in terms of 
the required resources such as number of ancilla levels D 
and possible operations on the atom-cavity system Vj^ . 

B. Arbitrary source-qubit interaction 

Now, we want to look at the problem from a more 
general perspective. We assume that the operators Vj^ 
are arbitrary isometries and that the ancilla decouples in 
the last step. We express the isometries in a given basis 

V=Y j V^\a,%m, (6) 

i,a>,(3 

where each V % is a D x D matrix and {|a), \/3)} are any 
of the D ancillary levels. This is the generalization of 
Eq. (j3]) and the isometry condition then reads 

l 

i=0 

The resulting n-qubit state is then given by 
l 

IV) = (<Pf\ V M--- V k\<Pi) Iw..,*i>. (8) 

ii ...i n =0 

This is a matrix-product state (MPS) [10, El| with D 
dimensional bonds and open boundary conditions, which 
are specified by the initial ancilla state \ipj) and the final 
ancilla state \<pf)- I n Figs. [2] and O we illustrate the 
sequential generation process schematically. 
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FIG. 3: Sequential generation of a multi-qubit state l^). In 
each step the D- dimensional ancilla produces one qubit. This 
process is described by a 2D x D dimensional isometry Vj^. 
The achievable multi-qubit states are instances of matrix- 
product states with D dimensional bonds and open boundary 
conditions specified by \cpi) and \<pf)- 



In practice, the question whether a given state can 
be generated with certain resources is more important. 
Therefore we have to show that every MPS of the form 



W) = ($f\%i]--'V[i]\<Pi), 



(9) 



with arbitrary maps V[k] ' Ha —> Ha ® Hb, can be gen- 
erated by isometries of the same dimension such that the 
ancilla decouples in the last step. Since the idea of the 
proof is an explicit construction of all involved isome- 
tries, it will provide a general recipe for the sequential 
generation of h^Y. Note that every state has a MPS rep- 
resentation [19l. [3Q[] . 
We start by writing 



({<p F \®l 2 )V [n] =V{ n] M b 



(10) 



where the 2x2 matrix Vj^ is the left unitary in the 
singular value decompositions (SVD) of the left hand side 
and M[ n ] is the remaining part. 

The recipe for constructing the isometries is the induc- 
tion 



(M, 



[fe] 



[fc-1] 



11]%!], (11) 



where the isometry V^^. x j is constructed from the SVD 

of the left hand side, and M[/ e _ 1 ] is always chosen to be 
the remaining part. 

After n applications of Eq. ([TTj) in Eq. (j9j), from left to 
right, we set \tpi) = M[i]| producing 



(12) 



Simple rank considerations show that Vj^_ fc j has dimen- 
sion 2 min [D, 2 k ] x min [D, 2 fc+1 ]. The dimension of the 
left unitary grows exponentially, i.e. V^_ fc+1 ] has di- 
mension 2 k x 2 fc , as long as 2 k < D. For 2 fc+1 > D, 
superfluous columns appear in V^ n _ k ^ since the original 

matrix (M[ n _ fc+1 ] 0I2) Vj n -fc] nas a ^ most D singular val- 
ues. Truncation leads to a 2 /c+1 x D dimensional isometry 

V [n-k}' 

Now M[ n _*.] has dimension D x D and all subsequent 
left isometries have dimension 2D x D. Therefore every 



Vjfcj can be embedded into an isometry Vjfc] of dimension 
2D x D. Physically, this means that we have redundant 
ancillary levels which we do not use. Finally, decoupling 
the ancilla in the last step is guaranteed by the fact that, 
after the application of Vj n _i], merely two levels of Ha 
are yet occupied, and can be mapped entirely onto the 
system Hb- This is precisely the action of Vj n ] through 
its embedded unitary Vj^ . 

Together with Eqs. (|5J) and (|9|), this proves the equiv- 
alence of three sets of n-qubit states: 

1. MPS with D-dimensional bonds and open bound- 
ary conditions. 

2. States which are generated sequentially and isomet- 
rically by a D-dimensional ancillary system which 
decouples in the last step. That is, the generation 
is deterministic. 

3. States which can be generated sequentially by a 
D-dimensional ancillary system in a probabilistic 
manner. That is, the preparation may only be suc- 
cessful with some probability and include measure- 
ments and conditional operations. 



In Sees. Ill CI and III D( we will show two other equivalent 
classes. We emphasize that this result holds as well for 
higher dimensional systems (beyond qubits) and that its 
constructive proof provides a recipe for the sequential 
generation of any state with minimal resources, namely a 
D-dimensional source/ancilla for a D-dimensional MPS. 

Note that in order to obtain MPS with periodic bound- 
ary conditions we require an interaction between the first 
and the last qubit. For systems, where an ancilla is 
available, one may store the first qubit within the an- 
cilla. Then, in order to produce MPS with D-dimensional 
bonds, we would require a 2D dimensional system to 
store the additional qubit. On the other hand many 
interesting states belong to the class of MPS with 2- 
dimensional bonds and in this case two additional atomic 
levels would suffice. 



C. The standard map 

The situation considered above assumes that arbitrary 
isometries can be achieved in the generation of a qubit, 
which amounts to have complete control over the source- 
qubit interaction. This does not quite correspond to what 
is the case in current cavity QED single-photon sources 
[4|, [5|, |9[ , where only atomic degrees of freedom can be ma- 
nipulated easily, while the isometry describing the gener- 
ation of a photon (qubit) is fixed. In the following we will 
show that also in this restricted scenario it is possible to 
generate arbitrary MPS with D-dimensional bonds if a 
2D-level atom is used as a source. 

We consider an atomic system with D states \ai) and D 
states as depicted in Fig.[H so that TLa — H a ®Hb — 
C D ®C 2 . 
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FIG. 4: Restricted interaction between ancilla and qubit. 
Each atomic transition from \ai) to its respective \bi) is ac- 
companied by the generation of a photon in a certain time-bin. 



That is, we will write for a superposition of 

\di) states, whereas |<p)|0) denotes a superposition of \bi) 
states. Since the last qubit marks whether the atomic 
level belongs to the \ai) or to the \b{) subspace, we will 
refer to it as the tag-qubit and write Ha — T~La' ® T~Lt- 

Now assume atomic transitions from each |a$) state to 
its respective \bi) state are accompanied by the genera- 
tion of a photon in a certain time-bin. This is described 
by a unitary evolution, from now on called "D-standard 
map" , of the form 



T: |^|1>t|0) s 
\p)M0)t\0) b 



\<p)a>\0)t\1)b , 
|^|0)r|0) s . 



(13) 
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FIG. 6: (a) Two-qubit gates V[{\ are sequentially applied be- 
tween nearest neighbor qubits. This situation can be sim- 
ulated by a two-dimensional ancilla as demonstrated in (b). 
Instead of applying the gate between neighboring qubits k 
and k + 1, we use qubit k and the ancilla and swap the ancilla 
state afterwards with the qubit k + 1. Since we can merge 
the swap operations and the arbitrary unitary operations V[{\ 
to arbitrary unitary operations, we know that the class of 
achievable states is equivalent to the class of MPS with 
open boundary conditions and two-dimensional bonds. 



Hence, T effectively interchanges the tag-qubit with the 
time-bin qubit. If, additionally, arbitrary atomic uni- 
taries Ua are allowed at any time, we can exploit the 
swap caused by T in order to generate the operation 

V\tp) = (0\tT(u a (\ V ) a ,\0)t)\0)b), (14) 

which is the most general isometry V : Ha' ~^ 7~La' ®T~Lb- 
Equation ([T4]) is also illustrated in Fig. [5l 



Ma 1 



\v)b ■ 

m-t- 



-|o)t 



> v\v)a< 



FIG. 5: The box represents a single step in the sequential 
generation scheme. Since tag- and time-bin qubit are always 
in state |0) after the preceding step, they can be ignored at 
the input. After an arbitrary unitary operation Ua on the 
atom, tag- and time-bin qubit are interchanged by the stan- 
dard map. Ignoring the tag qubit also at the output (since it 
is always in |0)), the process can be described by an arbitrary 
isometry V applied on the effective ancillary system A! . 



Therefore, the so generated n-qubit states include all 
possible states arising from subsequent applications of 
2D x D-dimensional isometries. On the other hand, they 
are a subset of the MPS in Eq. ([9|) with arbitrary 2D x D- 
dimensional maps, assuming that the atom decouples at 
the end. Hence, this set is again equivalent to the three 
mentioned above. 



D. Qubit-qubit interaction without ancilla 



In Fig. [6] (a) we depicted a system of N initially un- 
corrected qubits, which interact sequentially with their 
nearest neighbor qubit. 

This situation is in fact identical to the one consid- 
ered so far, as one can imagine the operation Vj^j being 
performed not between qubit k and k + 1 directly but 
between qubit k and a two-dimensional ancilla, which is 
then swapped on qubit k + 1 [see Fig. [6] (b)]. In the last 
step, the swap ensures that the ancilla decouples from 
the desired multi-qubit state Thus, also for direct 
qubit-qubit interaction, the class of achievable states 
is equivalent to the class of MPS with two-dimensional 
bonds (see [19| for a more formal argument). In section 
HV( we will come back to this scenario and consider an 
experimental setup, where always two neighbor atoms of 
a chain interact via a common cavity mode. 

If direct two-qubit interaction between neighbors is 
possible, there exists at least in principle - in contrast 
to the hen and egg scenario - no reason why one should 
not apply them more than once. If this is done m times, 
as illustrated in Fig. O still in a sequential manner, i.e. 
after an operation between the first and the last qubit 
one starts again with qubits 1 and 2, the achievable class 
of states will be described by MPS with bonds of dimen- 
sion d 2m_1 , that is for qubits 2 2m_1 . However, in this 
case the two sets will no longer be equivalent, i.e., there 
are MPS with D = 2 2m_1 for whose preparation we need 
more than m such layers of two-qubit gates. 
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FIG. 7: Repeating the procedure m (here m = 3) times 
defines a class of states described by MPS with 2 2m_1 - 
dimensional bonds. Note that in between the three sequences 
we also allow an arbitrary interaction between the first and 
the last qubit. 



FIG. 8: The D-\e\e\ atom can be viewed as a set of M qubits 
with D < 2 M . For an arbitrary operation we need local uni- 
taries for the atomic qubits and a universal two-qubit gate 
between the cavity qubit and one specific atomic qubit. 



III. SEQUENTIAL GENERATION OF 
PHOTONIC STATES 



There are two possible ways of encoding quantum 
information in photon pulses: orthogonal polarisation 
states and energy eigenstates. While the latter is more 
straightforward in terms of the required resources, po- 
larization encoding [3l|, [32[ avoids the trouble caused by 
a failure of the source for the generation of multi-qubit 
entangled states since a missing photon is not being mis- 
taken as an empty time-bin, i.e., the success of the en- 
coding step is heralded by the observation of a photon. 

The generation of multi-photon entangled states using 
several distant single photon sources was explored in Ref. 
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FIG. 9: (a) Atomic level structure: levels \a) (\b)) and |e) 
are coupled by a laser (cavity mode) off resonance, (b) After 
adiabatic elimination of the upper state |e), we are left with 
a Jaynes-Cummings type of Hamiltonian, where states \a,n) 
and |6, n + 1) are coupled. Both, the energy difference of 
those levels and the corresponding Rabi frequency depends 
on n. The reason for the first is the ac-Stark shift, whereas 
the second is due to the Jaynes-Cummings coupling. 



A. Time-bin entanglement 

We will first consider a situation where cavity decay is 
negligible on the time-scale of the operations performed 
on the atom-cavity system. After the cavity qubit has 
left, we start with the next step and repeating the pro- 
cess leads to a multi-qubit entangled photon state at the 
output as sketched in Fig. [TJ 



1. Arbitrary source-qubit operation 

In the following, we demonstrate how an arbitrary uni- 
tary operation on the 2D dimensional Hilbert space of 
the combined atom-cavity system can be realized. There- 
fore, we view the D-level atom as a set of M qubits with 
D < 2 M as depicted in Fig. [51 Then, we have to show 
that one can perform arbitrary two-qubit operations be- 
tween each pair of qubits, i.e., universal quantum com- 
puting. Since the atomic levels can be manipulated at 
will using Raman laser systems [34], [35|, it remains to 
propose an arbitrary two-qubit gate between one specific 
atomic qubit and the cavity qubit. 

Therefore, we consider a typical three-level lambda 
configuration (see Fig. [9J, where the hyperfine ground 
states | a) and \b) are coupled to the excited level \e) off- 
resonantly through a laser with Rabi frequency Q and 
detuning A + 5, and the cavity mode a with coupling 
strength g and detuning 5. 



Furthermore, we assume that the cavity decay rate k is 
smaller than any other frequency in the problem, so that 
we can ignore cavity damping during the atom-cavity ma- 
nipulations. In an appropriate interaction picture, the 
Hamiltonian of the system is then given by 



H = -A (a aa + a f a) + g(cr eb a + a) cr he ) 



~2^ e 



e m CTae) 



(15) 



with cr k i = \k){l\, {k,l} = {a, 6, e} and |A| > g, ft > 5. 
After adiabatically eliminating the excited state |e), the 
Hamiltonian for the effective D = 2 atomic system plus 
cavity mode, in an interaction picture with respect to 
-A(a aa + a 1 " a), is given by 
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(16) 

It describes an effective Jaynes-Cummings coupling be- 
tween the cavity mode and the atomic \a) — > |6) transi- 
tion with Rabi frequency gQ/2A. The other terms cor- 
respond to ac-Stark shifts. In an interaction picture with 
respect to the latter the Hamiltonian is given by 



ad 2^ 2A 

n=0 



c |a,n - l)(b,n\ + H.c. 

(17) 



Then, we choose the laser frequency such that 

-,2 



Q2 
4A 



A' 



(18) 
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For g 2 / A ^> gQ/2A, the effective interaction in all sub- 
spaces with n 7^ is then dispersive and the selective 
Hamiltonian [36[ is given by 

H sei = |^(|a,0)(M| + |M)M|) 

= + (19) 

where and erf denote the Pauli matrices acting on 
the atomic and the photonic qubit, respectively. Using 
a laser pulse of an appropriate duration, we obtain the 
entangling two-qubit gate 

VlSWAP = exp [i7r(|a,0)(M| + |6, l)(a, 0|)/4] . (20) 

Together with local operations in both qubits [37|, this 
suffices to generate an arbitrary two-qubit operation. 

2. Adiabatic passage 

In current CQED single-photon sources 0, [H, Ellij], an 
adiabatic passage is employed to realize the 1-standard 
map introduced above. Using one additional level, we 
will show how to generate familiar multi-qubit states like 
W [IK, GHZ [39L_and cluster states [40], which are all 
MPS with D = 2 [4l|. 

For this purpose, we consider an atom with three effec- 
tive levels {|a), |&i), I&2)} trapped inside an optical cavity. 
With the help of a laser beam, state \a) is mapped to the 
internal state |6i), and a photon is generated, whereas the 
other states remain unchanged. This physical process is 
described by the map 

M AB : \a) 1 * |&i)|l), 
|6i) - |&i>|0), 

\b 2 ) » \b 2 )\0), (21) 

and can be realized with the techniques used in Refs. 0, 
[1,0]. After the application of this process, an arbitrary 
operation is applied to the atom, which can be performed 
with Raman transitions. The photonic states that are 
generated after several applications are those MPS with 
isometries V[{\ = MabU^, with i = 1, . . . , n, being 
arbitrary unitary atomic operators. 

For example, to generate a W-type state of the form 

|^w) = e^ 1 sinei|0...01)+cos6ie^ 2 sin6 2 |0...010) 

+... + cos 61... cos B n -2e i ^" 1 sin 6 n _i |010...0) 
+ cos6i...cos6 n _i|10...0), (22) 

we choose the initial atomic state \cpj) = 1 62 ) and opera- 
tions U [ } = U h a \ 2 6i), with i = 1, . . . , n - 1, where 

U% ($i , 90 = cos &i\k) (k I + cos 9, 1 (I I + e 4 *' sin G t \k) (I \ 
-e- i ^sin6 i |/)(A:| + |ra)(ra|, (23) 



and {fc,Z,ra} = {a, 61,62}- To decouple the atom from 
the photon state, we choose the last atomic operation 
Ujfi = [/^ (0, 7r/2) and, after the last map M^g, the 
decoupled atom will be in state |6i). 

To produce a GHZ-type state in similar way, we choose 
\<p z ) = \a), = Ull^Q,), U [ } = Ull^/2), with 

i = 2, . . . , n - 1, and U% ] = Ufa (0, tt/2)^ (0, tt/2). 

For generating cluster states, we choose \cpi) = I62), 
U [ } = UX(*» ®i)U b a 2 bl (0, tt/2), with i = 1, . . . , n-1, and 
U [ a ] = KlS^&n)Ufa{^/2)U b a l^m, obtaining 

W=(g)(O?_ 1 |0) i + O?_i|l)i), (24) 

2=1 

where 0£_! = cos e»|0)i-i(0| - sin ©i|l>i_i<l| 

and 0|_! = e*** sin 8i|0)i_i(0| + cos © € |1>^_ 1 <1|, with 
z = 2, . . . , n — 1. Operators 0^_ x and act on the 

nearest neighbor-qubit i — 1 under the assumption 0[j = 
cos 61 and Oq = e z$1 sin 61. If one chooses <$>i = and 
Qi = 7r/4, this leads to the cluster states defined by 



B. Polarization entanglement 

Here, the cavity qubit is defined by single excitations in 
two modes a and 6, which have equal frequencies but or- 
thogonal polarizations. As above for the time-bin qubits, 
we will first show how to realize arbitrary source-qubit 
operations and then focus on the less demanding sce- 
nario, where the photonic qubits are generated by a stan- 
dard map, i.e., an adiabatic passage. 

1. Arbitrary source-qubit operation 

In this section, we will suggest how to implement an 
arbitrary operation on the atom-cavity system based on 
present cavity QED experiments 0, Therefore, 
one has to show how an arbitrary unitary operation on 
one specific atomic qubit {|a), |6)} and the cavity qubit 
{|l a ), (lb)} can be realized. We consider a double-lambda 
type atomic level configuration as illustrated in Fig. [lOl 
which couples to the two cavity modes via two Raman 
transitions. 

Two external laser fields drive the transitions from 
level I a') to the excited level \e a ) and from level |6') to 
the excited level |e&) with Rabi frequencies Vt a and f^, re- 
spectively. The cavity modes a and 6 couple to the tran- 
sitions between \e a ) and level \a) with coupling strength 
g a and |e&) and level |6) with coupling strength 

For large detunings 

|A a |,|A 6 | »0 a ,0 6 ,fi a ,ft 6 (26) 
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FIG. 10: Levels \a') and |6') (\a) and \b)) are coupled off- 
resonance to \e a ) and |e&) by two lasers (cavity modes). 



we can adiabatically eliminate the excited levels \e a ) and 
| eft) and end up with a Jaynes-Cummings type of Hamil- 
tonian, which is block separable in the subspaces spanned 
by {\a',n a ), |a,n a + l a )} and {\b',n h ), \b,n b + l b }} : where 
n a and n b denote the number of photons in mode a and 
b. As above, the ac-Stark shifts for n a = n b = can be 
compensated by choosing the frequency of the laser fields 
appropriately. Under the condition 



n a9l 9a^a ^bOj Qb^b 

A a >y 2A a ' A 6 ^ 2A b ' 
we obtain the selective Hamiltonian 



(27) 



-rj-p 9a^a i / n \ / -i | , 9b^b 



2A a 



2A 6 



|6',0 6 )(6,l 6 |+H.c., 



(28) 

where |0 a ) and |0b) denote the empty cavity modes a 
and b. For appropriate Q a (t) and f2f>(t), the evolution 
operator reads 

= exp [i7r(|a / ,0 a )(a, l a | + |6 7 , b ><6, U| + H.c.)/2]. 

(29) 

After the fcth photonic qubit leaked out of the cavity, 
the state of the system is given by 



\9 k ) = a\am)+0\bM), 



(30) 



where \if)%) and |^>£) are fc-qubit photonic states. Now 
we have to initialize the system for the next step, i.e., 
provide the polarization qubit 



-i(a\a',0 a )\r k )+P\b',0 b M)) 
a\a,l a )\r k )+f3\b,l b M), (31) 



where we applied U^ el in the second line. In order to real- 
ize an arbitra ry two-q ubit gate, we combine local opera- 
tions with a VIS WAP two-qubit gate between the atomic 
and the photonic qubit. The latter can be achieved in 
three steps, 

VlSWAP = (U^y 1 j*(\"'W\+\b')(a'\)/* v ^ (32) 

It remains to show how to decouple the atom from 
the generated multi-photon state in the final step n. 
Therefore, one has to map the atomic state on the last 
photonic qubit: after the transformation of the atomic 
levels | a) — > — i\a') and \b) — > —i\b'), driving only the 




FIG. 11: Atomic level structure: level \a%) (\bi)) are coupled 
to \fi) via an adiabatic passage driven by a laser with Rabi 
frequency Q a (£l b ) and the cavity mode a (b) with coupling 
strength g a (g b ). 



\a>')\0) — > |«)|la) transition by choosing Q b = in U^ el 
leads to 



|*n- 



a\a,l a )\i> a )-i(3\b',0 b )\i> b ), (33) 



where the cavity mode b remains empty and level \b) is 
not yet populated. Now we transform \a) — > |6), then 
apply U^ e i again and, since the photon in the mode a can 
not be absorbed, end up with 

|* n > = |6>®(a|l >|^ )+ ) 9|l 6 >|^)) ) (34) 

where the atom decouples in its final state \b). 

2. Adiabatic passage 

We consider an effective three-level system, with 
atomic ground states |a^), \bi) and |/^), which couples 
to the two cavity-modes a and b through two indepen- 
dent adiabatic passages controlled by two external lasers 
with Rabi frequencies Q a and Q b as depicted in Fig. [TTJ 

In each generation step the standard map is achieved, 



\ a i) 



\fi)\U) 

\mu) 



\h)\la), 



(35) 



where we first applied the adiabatic passages and then a 
unitary operation \fi) — > \bi). In general, we have i = 
1,...,D. Instead of an effective 2D-\eve\ atomic system, 
as for the generation of time-bin entangled MPSd, we 
require a 3D-level atomic system. 

Note that the results for the generation of time-bin 
qubits apply: for W-type, GHZ and cluster states we 
require only D = 1 and one additional level, i.e., an ef- 
fective four-level atom. 



IV. SEQUENTIAL GENERATION OF ATOMIC 
STATES 

For the generation of atomic multi-qubit states the cav- 
ity mode and the atom interchange their function. Now, 
the cavity field is employed as the ancillary system to 
interact with initially uncorrelated atoms which sequen- 
tially pass through the cavity. Therefore, we require a 
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FIG. 12: A stream of uncorrelated atoms crosses a cavity. The 
atomic qubits couple sequentially to the cavity mode, which 
acts as a 2-dimensional ancillary system. If an arbitrary uni- 
tary operation can be realized between the cavity qubit and 
the atomic qubits, the class of entangled multi-qubit atomic 
states at the output are equivalent to the class of MPS2. 



very stable cavity field which is provided typically by 
microwave cavities [28|, [29|. In the second part of this 
section, we consider another scenario where the atomic 
qubits interact directly via a common cavity mode [421 ]. 



A. Cavity acts as ancillary system 

In the first scenario the atoms pass through the cavity 
in such a manner that only one atom couples to a single 
cavity mode at a time, as depicted in Fig. [121 We employ 
the same atomic level configuration as in section 3A and 
define the ancilla qubit as the cavity mode Fock states 
|0) and |1). The atom-cavity system is then described 
by the selective Hamiltonian from Eq. [19] Choosing the 
field frequency and pulse length appropriately leads to an 
IS WAP gate between the photonic and the atomic qubit. 
In the basis {|6, 0), |6, 1), \a, 0), \a, 1)}, 



ISWAP 



( 1 

i 

i 

\ 1 . 



(36) 



Three applications of the ISWAP gate accompanied by 
local unitary operations suffice for an arbitrary two-qubit 
gate [27], though they could only be applied in setups 
where the position of the atoms is fully controlled [y, [j], 
Q. In fact, it would not only allow the generation of 
all MPS with two-dimensional bonds, but also all other 
states since the atoms could be moved back and forth, 
allowing for universal quantum computing. One of those 
setups use optical cavities and, unfortunately, the cavity 
decay is not negligible on the time-scale in which the 
atoms are moved in and out of the cavity. 

For setups with high- finesse microwave cavities [28], [29| , 
the atoms pass once the cavity and local rotations can be 
implemented via suitably placed Ramsey zones. So the 
natural question arises: which states can be generated if 
only one ISWAP gate and local unitaries on the atomic 
qubits are available? 

If all atoms are prepared in state \a) and the cav- 
ity mode is initially empty, sequential application of the 



VlSWAP gate leads to 



i|0)|a 2 ,ai) + ^|l)|62,ai) + -^|l)|a2,6i) 



|a n _i . . . ai) 



+ ^-_ |l)|& n _i,a n -2, . . . ,ai) + . . . 
i 

+ -y=\l)\a n -i, . . . ,a 2 , 61). 



(37) 



This W-type state is still entangled with the cavity qubit. 
One way of solving this problem is to measure the cavity 
qubit in the basis (|0) ± \ l))/>/2. 

Let us concentrate on the generation of the one- 
dimensional cluster state given in Eq. [25j The required 
control Z (CZ) gate between neighboring qubits can be 
realized through a CZ gate between the ancilla and the 
first qubit followed by a SWAP operation on the ancilla 
and the second qubit. Therefore, we decompose 



ISWAP = 



/ 1 

10 

10 

V 1 



/ 1 

i 00 

i 

\ 1 



(38) 



where the first matrix represents the SWAP operation 
and the second matrix is equivalent to the CZ gate up 
to local operations. As has been depicted in Fig. [13] (a), 
the ISWAP can then be written as 



ISWAP = i SWAP CZ [R z (ir/2) <8> R z (tt/2)] , (39) 



with 



R z {4>) 



e^/ 2 j ■ 



(40) 



Since it is difficult to perform local operations on the 
cavity qubit, say in step fc, we make use of the SWAP 
and apply it instead in step k + 1 after the SWAP on 
qubit k + 1. Fortunately, R z (tt/2) commutes with the 
CZ gate and we obtain the recipe shown in Fig. [13] (b). 
The last unitary on the ancilla can be ignored because 
it has no influence on the desired multi-qubit state. In 
order to obtain the cluster state as defined in Eq. [25l we 
have to compensate for the local unitaries on the atoms 
and prepare all atoms initially in superpositions of (|0) + 

|1»/V2. 



B. No ancilla: atoms interact via cavity mode 

Using the cavity mode as an ancilla has the disadvan- 
tage that it may decay during the time-interval between 
two successive atoms. This problem can be avoided by 
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(a) 



(b) 



^U-t-j- 



1S-m- 



FIG. 13: Cluster state generation, a) An ISWAP gate is per- 
formed sequentially between the atomic qubits and the cavity 
mode, b) One can as well assume that the local unitaries 
are performed on the atomic qubit k after the fcth SWAP. 
Therefore they can be compensated in a final step. 




Product state 



FIG. 14: A stream of uncorrelated atoms passes a cavity. Two 
atomic qubits couple to each other via the common cavity 
mode. If an arbitrary unitary operation can be realized be- 
tween them, the class of entangled multi-qubit atomic states 
at the output are quivalent to the class of MPS 2. 



a direct interaction between the atomic qubits via the 
common cavity mode. As sketched in Fig. [HJ there 
are always two atoms at the same time inside the cav- 
ity. Zheng and Guo [43[ proposed the implementation of 
a vISWAP between the atomic qubits via a dispersive 
scheme. Note that in order to implement a CNOT an 
additional atomic level would be required [43[. Recent 
experiments [7|, [9| raise hope that it will soon be possible 
to move two atoms into a cavity in a well controlled man- 
ner. Then the vISWAP accompanied by local unitaries 
suffices to perform an arbitray two-qubit operation. 

If we consider that both atoms cross the cavity only 
once [44j , where an vISWAP has been demonstrated, the 
question arises again: which mutli-qubit atomic states 
can be sequentially generated in this manner? Since the 
atomic qubits can be manipulated locally before and after 
the gate with Ramsey zones, we gain more possibilities 
than in the corresponding ancilla case. 

The generation of the W and cluster states follows the 
lines of the previous case of secion IVB, since we only 
replace the ancilla by the neighbour qubit and the su- 
perfluous SWAP gate does not affect the desired output 
state. In order to engineer a GHZ state, we use another 
decomposition of ISWAP gate, given by 



ISWAP = zSWAP [R z (7T/2)^R z (ir/2)} 
x [1 <8> H] CNOT [1 <8> H] , 



35: 



FIG. 15: The ISWAP gate between neighboring qubits is 
equivalent to a SWAP followed by a CNOT up to local uni- 
taries. The latter can be compensated since local operations 
can be applied to the atomic qubits with Ramsey zones. 



where the Hadamard gate is given by 



y/2 



1 1 
1 -1 



(42) 



Since we can compensate the local unitaries acting on the 
atomic qubits, we end up with a sequential application of 
a CNOT followed by a SWAP on all neighboring qubits 
as depicted in Fig. [151 

Now, let us assume that our initial state has all 
atoms prepared in level \a), except the first which enters 
the cavity in the superposition (\a) + \b))/y/2, 



T2 



\a,...,a) + \b,...,b)). 



(43) 



(41) 



Here, the arrows in step i indicate the application of a 
CNOT and a SWAP between qubits i and z+1, producing 
at the end the desired GHZ state. 

Further experimental efforts will make use of two con- 
secutive cavities in the same setup [45|. This opens up 
new possibilities in terms of state engineering in the light 
of the present work. Three cavities in a row, which are 
crossed by atoms such that unitary operations may be 
performed between the cavities, would even allow for 
arbitrary two-qubit operations between the atoms and 
therefore lead to a class of states equivalent to the class 
of MPS 2 . 



V. CONCLUSIONS 

We developed a formalism to describe the sequential 
generation of entangled multi-qubit states. It classifies 
the states achievable by a D-dimensional ancilla as MPS 
with D-dimensional bonds and provides a recipe for the 
sequential generation of any state. It turns out that all 
states that can be generated non-deterministically with 
a D-dimensional ancilla also belong to the class of MPS 
with D-dimensional bonds. Therefore, we will be able to 
provide a recipe for their deterministic generation as well. 
Remark that the formalism applies also to a situation 
where qubits interact directly in a sequential manner. 

For the generation of entangled multi-qubit photonic 
states by a cavity QED single-photon source, we propose 
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am implementation of an arbitrary interaction between 
the atom and a cavity qubit defined either by the absence 
and the presence of a photon or by the polarization of 
single excitation in the cavity. Moreover, we discussed 
the case of an adiabatic passage for the photon generation 
as being employed for single-photon generation in current 
cavity QED experiments 0, [H, 0] . 

For coherent microwave cavity QED experiments [Hj], 
where atoms sequentially cross a cavity and interact with 
the same cavity mode, we give a recipe for the generation 
of W-type states as well as cluster states. We considered 
also the case of direct coupling of successive atoms via 
the cavity mode. For this case we show how to generate 
GHZ states, in addition to the W-type and cluster states. 

Other physical scenarios as a light pulse crossing sev- 
eral atomic ensembles [46| or trapped ion experiments, 
where each ion interacts sequentially with a collective 
mode of the motion |33, [35|, [47j , may also be described 



by the present formalism. 
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